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Abstract: We decompose the quark propagator in the presence of an arbitrary gluon 
field with respect to a set of Dirac matrices. The four- dimensional integrals which arise 
in first order perturbation theory are rewritten as line-integrals along certain field lines, 
together with a weighted integration over the various field lines. It is then easy to trans- 
form the propagator into a form involving path ordered exponentials. The resulting 
expression is non-perturbative and has the correct behavior under Lorentz transforma- 
tions, gauge transformations and charge conjugation. Furthermore it coincides with the 
exact propagator in first order of the coupling g. No expansion with respect to the in- 
verse quark mass is involved, the expression can even be used for vanishing mass. For 
large mass the field lines concentrate near the straight line connection and simple results 
can be obtained immediately. 
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1 The quark propagator 



The quark propagator S{x, y; A) for a quark of mass m in the presence of a gluon 
field A'^ plays an important role in many investigations of quantum chromodynamics. 
It appears, e.g., if one considers the quark four-point Green function, the basis of all 
modern investigations on quark-antiquark interactions, and integrates over the quark 
fields. It is defined by Q 

S^f3{x,y-A) = -i< 0|T(V^,(x)V^^(i/))|0 > . (1.1) 

The spinor indices a, (3 will be dropped in the following. We recall the relevant properties 

of the propagator. With the covariant derivative D^{x) which acts on operators to the 

< — * 

right, and (y) which acts on operators to the left, defined by 

d 

Kiy) = ^ + '9My), (1-2) 

the field equations give 

[tYD^{x)-m]S{x,y;A) = 6^'\x-y), 
Six,y;A)[-tYDliy)-m] = 5^'\x-y). (1.3) 
These equations may be reformulated as integral equations: 



S{x,y]A) = So{x-y)- g J S{x,z;A)'y''Af,{z)Soiz-y)d'^z 

= Soix -y)~gj S,{x - z)^^ A,{z)S{z, y- A)d\ (1.4) 
with 5*0 the free propagator. From charge conjugation one has 

S{x, y- A) = 7'7°5^(y, x; -A^hV- (1-5) 
Finally, under a gauge transformation ip ^ ip' = e''®ip,ijj Tp' = 'ijje~^^,A^ A'^ = 
e*®74^e~*® — (2/(7)((9^e*®)e~*®, the propagator transforms as 

S{x, y; A) S'{x, y; A) = e*®(^)5(x, y; A)e-'^^y\ (1.6) 

An exact solution for S{x, y; A) for an arbitrary gluon field A^ is not available. However, 
it would be highly desirable to have an approximation which respects the fundamental 
properties of the propagator, in particular the correct transformation under Lorentz 
transformations and under gauge transformations. 

The two well known approximations, perurbation theory and the static approxima- 
tion, either violate gauge covariance or Lorentz covariance: 

^Our formulae are given in Minkowsi space, we use Bjorken Drell conventions Q, and the field tensor 
is defined by F^i, = df_,Ai, - d^Af, - ig[A^, A^]. 
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Perturbation theory: Iteration of eq. (1.4) gives the perturbation series 

S{x, y; A) = So{x - y) - g J So{x - z)YA,{z)So{z - y)d'z + ■■■. (1.7) 

Any finite order of the perturbation series gives the correct behavior under Lorentz 
transformations and under charge conjugation. It will, however, never be able to 
describe central features of QCD like confinement. In particular, it clearly violates 
gauge invariance; no truncation of the perturbation series (1.7) has the correct 
transformation property (1.6). 

Static approximation: Following the pioneering work of Brown, Weisberger p| and 
Eichten, Feinberg [Q, one neglects the spatial part i'y"^Dm{x) in eq. (1.3). The 
partial differential equation then becomes an ordinary differential equation which 
can be solved in closed form. This leads to the static approximation 

r 1 + ^0 ^ _ ^0 

Sstat{x,y;A) = ie{xo - yo) — ^ <S>{yo - Xq) — - — 

(5(=^)(x-y)e-''"i"o-^«IPexp{i^ f A^{z)dz'^}. (1.8) 

The path in the line integral is the straight line from y to x, and the path ordering 
orders Aq{x) to the left, ■ ■ ■ , Ao(?/) to the right. 

The neglected spatial term i^j^^Dmix) in (1.3) can subsequently be taken into 
account as perturbation. This approach has been extremely succesful (see e.g. 
the review [^). It can be easily generalized to quarks moving with four velocity 
v^, and thus is the direct predecessor of heavy quark effective theory. Succesful 
combination with perturbation theory has also been made more recently (see 
and references therein). For quarks moving with high momentum a related formula 
can be derived from the eikonal approximation 0]. 

The static approximation and it's generalizations are non-perturbative and have 
the correct behavior under gauge transformations and under charge conjugation. 
However, they drastically violate Lorentz invariance. Therefore the static approx- 
imation is useful for heavy quarks only. It needs quite an effort to recover the 
relations following from the original Lorentz invariance subsequently |Q, 0] . Finally 
the static approximation does not coincide with the exact propagator even in the 
trivial case of vanishing gluon field. 

It is obvious that a gauge covariant propagator should contain path ordered exponen- 
tials as in (1.8). Using just the path along the straight line between x and y would give a 
Lorentz covariant result, but such a procedure would be far too simple. It would involve 
the vector potential only along the straight line connection and nowhere else, which is 
clearly unphysical. We prefer to proceed systematically by rewriting perturbation theory 
in a suitable way. It can then easily be transformed into a gauge covariant expression 
by simple exponentiation, while keeping the correctness of perturbation theory in the 
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relevant order. The non-perturbative approximation for the quark propagator obtained 
in this way has the following properties: 

• Correct behavior under Lorentz transformations. 

• Correct behavior under gauge transformations. 

• Correct behavior under charge conjugation. 

• Agreement with perturbation theory in first order of the coupling. 

The representation is a weighted superposition of path ordered exponentials between x 
and y along well defined field hues. We don't need to assume that the quark mass is 
large, we could even put it equal to zero. This opens perspectives to applications which 
were hard to attack previously. 

The paper is organized as follows: 

In sect. 2 we decompose the propagator with respect to Dirac matrices, and write 
the formula of first order perturbation theory in a form which is convenient for the fol- 
lowing. In sect. 3 we rewrite the four-dimensional space-time integrals which arise in 
perturbation theory as weighted superpositions of line integrals over certain field lines 
which all run from x to y. From this form one can simply derive a representation in 
terms of superpositions of path ordered exponentials. This representation coincides with 
perturbation theory up to order g and has the correct behavior under gauge transfor- 
mations. In sect. 4 we evaluate the weight function explicitly and show a plot of the 
field lines for different masses. We discuss the limit of large mass m, and give some first 
simple applications. Actual applications will be given in forthcoming papers. 

2 A useful form of first order perturbation theory 

We start with the first order approximation (1.7), and express the free propagators 5*0 
by the free scalar propagator A in the following way: 

Sq{x - z) = A{x - z)[-iji, d /dzi, + m], 

So{z-y) = [t-fxd/dz^ + m]A{z-y). (2.1) 

The free scalar propagator satisfies the equation 

{d,df + m^)A{x) = -5(^)(x). (2.2) 

In the following all differential operators in the integrand are understood as differentia- 
tions with respect to the variable z. 

Using the well known identities 7^7^ = g^f, - ia^,^, and 7j.7^7a = g,,i^^\ - guxltj, + 
QixXlv + i^vuXnl^l'^ 1 one can write the propagator in form of the familiar decomposition 

S{x, y;A)=s+ rf' + v'^'y^ + a^'^y^ + t^'a^,. (2.3) 
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Before giving the expressions for s, ■ ■ ■ , t^'^ which result in this way, it is convenient for 
later use to define a scalar field u{z;x,y) and a vector field u^{z;x,y) by 

u{z;x,y)^A{x-z)A{z-y), (2.4) 
and (with d /dz^- d /dz^) 

ui'iz; X, y) = - A-^ {x - y) [A(x - z) ^ A(z - y)] . (2.5) 
Prom (1.7), (2.1), and (2.3) we then obtain the following equations: 

s = mA{x-y)[l + ig J u^'{z;x,y)A^{z)d''z], (2.6) 
p = 0, (2.7) 
= '-d/dx^A{x-y)-'-A{x-y)d/dyi, (2.8) 

-g J |a(x - z)[m^g''''+ - q"" f + d''t]A{z - y)} A,{z)d^z, 

^ igei^-^- J[A{x-z)dJxA{z-y)]A,{z)d''z, (2.9) 
J {dMz; X, y)A,{z) -ifi^ u)} d'z. (2.10) 



We next transform in (2.8), we denote the four terms in the second line by Vi + 

+ "^3 + "^4- "^2 can use dz= ~ dy, in correspondingly dz= ~ d^- The 
differentiations with respect to y and x can then be taken outside of the integral. The 
term is split into two identical contributions, in the first one we use eq. (2.2) for 
A{x — z), in the second one for A{z — y). Next perform a partial integration on one of 
the derivatives in the d'Alembert operator. The terms where the differentiations act on 
the other A cancel against the contribution i>3 , and one has the intermediate result 

v^ + v^ = ^A{x-y){A^ix)+A^iy)} 

-| J dxu{z;x,y)[d^A>'{z) - d''A\z)]d^z 

-| J dxu{z; X, y)d^A\z))d'z. (2.11) 

We subtraced and added the term d^A^{z). In the second line we can then replace 
d^A'^{z) — d'^A^{z) by F^^{z)^ which only introduces a higher order error 0{g^). Fur- 
thermore the differential operator d\ = d\ which acts on u{z;x,y) can be replaced by 
—{d^ + dy) and taken outside of the integral. In the third line we perform a partial 
integration on 9^, shift the differentiations from the variable z to x and y, and take the 
differentiations out of the integral. After these manipulations can be expressed in 
a rather compact form if wc replace partial derivatives by covariant derivatives which 
introduces a correction of order g^ only: 
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= '-D^{x)\^Aix-y)[l + igJu''iz;x,y)A,iz)dh]^ 

|a(x -y)[l + ig Ju''{z;x,y)A,{z)d'z]^ D*' (y) 
-|l>,(x) J u{z;x,y)F^^{z)d'z - |/ u{z;x,y)F^''iz)d'z dI (y). (2.12) 

In we perform a partial integration with respect to 9k or 9a ; take half of the sum of 
both terms, antisymmetrize the dKA^{z) or dxA^{z) in the integrand, and introduce the 
field strength tensor F^u{z) or F\y[z) as before. Shifting again the differentiation from 
^ to x and y we have in order g 

= J e'^""^ [Dx{x) j u{z-x,y)FUz)dS - j u{z; x,y)F,^{z)d^z D, (y)] . (2.13) 
Similar manipulations applied to t^'^ finally give 

^ y u{z; X, y)F''''{z)d''z. (2.14) 

The decomposition (2.3), together with the formulae (2.6), (2.7), (2.12), (2.13), (2.14) 
is now in a form which allows to rewrite the four-dimensional integrations d'^z as a 
superposition of line integrals. 

3 Gauge covariant reformulation 

The vector field u^{z;x,y) defined in (2.5) satisfies the fundamental equation 

^}^!li^l^^S^'\z-y)-S^'\z-x), (3.1) 

which is a simple consequence of eq. (2.2). Therefore may be interpreted as a four- 
dimensional velocity field of an incompressible fluid with a point-like source at y and 
a sink at x. The stream lines z'^{s,w), which all run from y to x, are defined by the 
characteristic equations 

'-^^^WUM). (3.2) 

Here s is the parameter which describes the motion along the stream line, while the three 
dimensional parameter set w characterizes the various stream lines. To make s unique, 
it is convenient to fix s = at the symmetrical point of the stream line which has equal 
distance to x and y. The dependence on x, y has been suppressed in the notation. There 
is precisely one field line passing through every space time point, except, of course, for 
the source points x and y. After having solved (3.2) there is a unique correspondence 
between the 4-dimensdional space-time coordinates z'^ and the parameters s, w, i.e. 
(2;°, z^, z^, z^) ^ (s, w^, w^). 
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We next write the four-dimensional integrals over d'^z which appear in s and the first 
two hnes of (eqs (2.6) and (2.12)) as integrals over dscPw, the Jacobian is called p(w): 

o[s, w'-, w^, w^) 

We have anticipated the crucial fact that p does not depend on the curve parameter 
s. This is a direct consequence of the incompressibility of the flow, and easily proved 
from the following geometrical argument. Take an infinitesimal four-dimensional box 
in (s, w)-space with corners {sq, w) and (sq + ds,w + dw). In z-space this corresponds 
to a an infinitesimal region with a certain volume. Consider now the motion of this 
volume along a field line from sq to Si, keeping w,dw and ds constant. Because of the 
vanishing divergence (3.1) outside the sources, the volume in 2;-space stays constant, the 
volume in (s, w)-space stays constant anyhow by construction. This demonstrates that 
the Jacobian (3.3) is indeed independent of s. 
We can therefore write 



lu'^(z;x,y)A,(z)dh = J p(w)[l^^ A,(z(s,w))u'^(z(s,w))ds]d' 

= J p{w)[£Adz]d^w. (3.4) 

RereyfyAdz is a shorthand notation for the line integral of from y to x along the 
stream fine characterized by the parameter w. 

The normalization of p(w) is easily obtained from the special case A^j^^z) — dij,Q{z), 
where both sides of (3.4) can be immediately integrated. This gives 

J p{w)d^w = 1. (3.5) 
We are now in the position to rewrite the scalar function s in (2.6): 



s — mA{x — y) j p{w)[l + ig Adz]d^w (3.6. a) 

= mA[x — y) J p{w)P exTp[ig ^ Adz]d'^w + 0{g'^). (3.6. b) 

This was the essential step of the approach! Once having written the first order ap- 
proximation in terms of line integrals, exponentiation allows to promote it into a non- 
perturbative gauge covariant expression. 

One may wonder about the justification of this step which led to a non-perturbative 
expression simply by exponentiation. It does not make sense to compare the two lines 
in (3.6), because (3.6.b) is gauge covariant while (3. 6. a) is not. In fact one can easily 
see, e.g. for a suitable pure gauge, that (3. 6. a) can be made as large as one likes. 
The correct question to be asked is whether (3.6.b) is a reasonable approximation to 
the exact propagator. To answer this question one only needs to check the quality of 
the approximation in some special convenient gauge. The manifest covariance of the 
expression then guarantees this quality for any gauge. A sufficient condition is the 
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smallness of g^yAdz for all relevant stream lines, which justifies the validity of first order 
perturbation theory. Clearly one can always find a gauge where Af^dz'^ = on one special 
stream line, but it is not possible to do the same for two or more lines. For heavy quarks 
the situation is simple. We will show in sect. 4 that in this case all relevant paths lie 
near the straight line connection. If one chooses a gauge such that A^{x — y)^ vanishes 
along this line, and if the variations of A^ near this line are small, the approximation will 
be justified. In applications we have to perform an integration over the gauge field A^ at 
the end, which involves the whole spectrum in momentum space. Because there are only 
two scales available, m and Aqcd, one expects a good approximation if m >> Aqcd- 

For light quarks the lines spread out over the whole space, and the above argument 
cannot be applied. But in any case the step of exponentiation leading from (3. 6. a) 
to (3.6.b) can be interpreted as partially including higher orders of the perturbation 
series, namely a minimum of those necessary to guarantee gauge covariance. There are 
good reasons to believe that gauge covariance is such a fundamental principle that it 
may indeed be used to transform perturbative expressions into non-perturbative ones of 
physical relevance. 

The expressions in the first and second line of f ^ can be treated in exactly the same 
way. The integrals in the third line of f^, as well as those in and t^'^ have a different 
form, but the structure of all of them is identical. They have a factor g in front, therefore 
one can introduce path ordered exponentials without changing the result in order g. To 
each z ^ x,y there belongs a unique s' and w which characterize it's position s' on the 
field line w. One can write 

gF'^^izis',w)) = P !^gF'^^izis',w))exp[tg£Adz]'j + Oig'). (3.7) 

We have chosen the symbol s' in z{s',w) in order not to mix it up with the curve 
parameter s in the path ordered exponential. The color matrix F'^'^{z{s' has to be 
included in the path ordering prescription with respect to s of the field line characterized 
by w. In this way (3.7) behaves correctly under gauge transformations. 

Finally one has to multiply by u{z;x,y) and perform the integration over d'^z = 
p(w)ds'd^w, resulting in 

J u{z-x,y)gF'^''{z)dH = 

j p(w)u(z(s',w))p|^F^'^(z(s',w))exp[«^^Ad2]|ds'd3'u; + 0(/). (3.8) 

This type of integral is a generalization of the operator insertions into a Wilson loop 
introduced by Eichten and Feinberg and later on widely used in the hterature. It is, 
however, more general, because the insertions are not equally distributed along the path, 
but weighted by the s'-dependence of m(z(s', w)). 

We summarize our representation for the quark propagator: It has the decomposition 
(2.3) with 
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s — mA(x — y) J p(w)P exp[ig -j/ Adz]d^w, (3.9) 
p = 0, (3.10) 

'f" = ^D''{x)I^A{x-y) J p{w)Pexp[ig^''Adz]d^w^ 

|a(x -y)l piw)Pexp[ig£Adz]d^w^ D*'' (y) (3.11) 
~D^{x) J p{w)u{z{s' ,w))P l^gF^^izis' ,w)) ex.p[ig ^^Adz]^ ds'd^w 
~ J p{wHz{s',w))P l^gF^%z{s',w))exp[ig£Adz]j ds'd^w dI {y), 
^e^^^'^{L',(x) J p{w)u{z{s\w))pf^gF,Ms':^))exp[igfAdz]^ds'd'w 



- J p{w)u{z{s',w))P^^gF,Ms',w))exp[ig£Adz]^ds'd^w d\ (?/)}, (3.12) 
= j p{v^)u{z{s\w))P^gF^'^{z{s\w))Pcxp[tg£Adz\^ds'd'^w. (3.13) 

We could have simplified the curly brackets {• • •} in v*^ by s/m, but we left it in the 
present form in order to keep it applicable also in the case of small or vanishing mass. 
Obviously the last two lines in and the terms and t^^ contain the same types of 
insertions into path ordered integrals. 

Clearly the representation (2.3), (3.9)-(3.13) fulfills all properties mentioned at the 
end of sect. 1. 



4 Weight function, stream lines, and a simple appli- 
cation 

It is now appropriate to switch to Euclidean space, i.e. put = —ixf, x^ = x^, 7° = 
74^! 7" = ^7rf • III the following the index E will be written explicitly only where it appears 
appropriate. For the explicit calculations we make use of the rotation symmetry around 
the vector {x — y)^. Choose, just for intermediate simplification of notation, a system 
where {x — y)^ has a four-component only and where x = y = 0. It is then convenient to 
describe the vector 2;^ by it's four-component 2:4, and ordinary three-dimensional polar 
coordinates r, ©, </?, i. e. 

2;^ = (r sin cos (/?, r sin sin y?, r cos 0, Z4). (4-1) 

Let us now choose an appropriate parametrization of the various stream lines. We clas- 
sify them by the orthogonal distance w of the line from the midpoint [x + y) /2 between 
the sources, and by the angles and </?, thus d^w — dwdQd(p. Fig. 1 shows a stream 
line together with the parameters introduced above. 
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Fig. 1: The plane spanned by the stream line, and the coordinates introduced in the 
text (in the system where x — y has a four- component only). 

The weight function p(w) defined in (3.3) becomes 

^i^' ^) - e, z,) d{s, w, d{s, w) ■ ^^-^^ 

It is now very convenient that p depends on w only, but not on the curve parameter s. 
Therefore we can evaluate it at a suitable point. We choose the symmetry point ( in the 
middle of the stream line with the coordinates r = w , Z4 = (x + y)^/'^-. At this point we 
obviously have dr/ds = and dr/dw = 1, and thus d{r, Z4)/d{s,w) — —dz^^/ds = —U4, 
where we used the definition of the stream lines in (3.2). This allows to write down the 
weight function in closed form. 

p{w,Q) — —w'^ sinQ U4^{(;x,y), (4.3) 

with 

Cm = ("^ sin cos w sin 6 sirup, w cos 0, {x + y)i/2). (4.4) 
Obviously u^^C,; x, y) is independent of the angles 6, ip. 



It is convenient to put 



p{w,e) = P^smQ, (4.5) 



such that the weight function p{w) is normalized to 

roo 

/ p{w)dw = 1. (4.6) 
Jo 

We first give the resulting formulae for the special case of vanishing quark mass which 
may be of some general interest: 

(4.7) 
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{z, X, y) 27^2 I _ ^)4(^ _ + _ ^)2(^ _ 1 ' l^.yj 

p(°)H^ ^"^''f ,.3 - (4.10) 

The maximum of p^^^ is at w!^\^ = \x — y\/{2\/2). 

We now come to the general massive case. The free scalar propagator A then is 

A(x) = -^K^{mx), (4.11) 

which gives 

u{z; X, y) = Jll_M^\^ --\)Ki{m\z-y\)_ ^^_^2) 



With the relation {Ki{z)/z)' = —K^iz^jz for the Kelvin function, one further obtains 



2 

m\ 


\x - 


-y\ 




i K2{m\ 


\x — z\)Ki{m\ 






[m\ 


\x - 


-y\)\ 


[ {x — zy\ 


\z-y 





\z-y)\. (4.13) 



Kiim\ 


a; — 2; 


)K2{m\ 


\z-y\) 




x — z\ 


\{z-yY 



m^w\x - yf K,{msJ{x - y) 74 + w^)K2{msJ{x - y)V4 + w'') 
^^'^^ TT Xi(m|x-y|)[(x-y)V4 + «;2]3/2 " ^^'^^^ 

The weight function p(w) is trivially suppressed for small w by the volume element in 
polar coordinates, it rises to a maximum at some Wmax, and decreases exponentially for 
large w. 

In Fig. 2 we plot the function p{w) for fixed distance |a; — 1/| for various values of the 
mass. It is clearly seen how the maximum moves to the left if the quark mass increases. 
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w/\x - y 



0.2 0.4 0.6 0.8 1 1.2 1.4 

Fig. 2: The weight function p{w) defined in (4-6). The distance \x — y\ is fixed, 
we show the curves for the mass values m\x — y\ — 0, 1, 5, 10, 20. The maximum moves 
from right to left and increases with increasing mass. 

In Fig. 3 we show the stream hnes of the vector field for four values of m\x — y\. 




-0.4 - -0.4 - 

Fig. 3: The stream lines (3.2) for the cases (from left to right, top to bottom) 
a) m\x — y\ = 0, h) m\x — y\ = 5, c) m\x — y\ = 20, d) m\x — y\ = 100. The 
sources at x and y are located at ±0.5. We show the stream lines for the values of 
W/Wmax = 0, ±0.5, ±1, ±1.5. 

For vanishing mass the lines spread out widely in space, up to w of the order of\x—y\. 
For increasing mass they concentrate more and more to the straight line connection. 
Apparently the product m\x — y\ has to become quite large, however, in order to get a 
sizeable concentration. 

The large mass limit will now be investigated analytically. We use the asymptotic 
behavior Kn[z) ~ a/vt/ {2z)e~^ for 2; — >• oo and get 
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^^^^^4(2^^^""^' form^oo, (4.15) 

^ ^-m[\x-z\+\z-y\] 

u{z;x,y) —^j — form^oo, (4.16) 

SzTT"* [\x — z\\z — y\\'^'^ 



3/2 

o \ o_i„ -,11-, „.i ; I |x — \z — y 



m\ 


\x - 


-y\ 




27r| 
m 


\x - 


- z\ 
oo, 


\z ■ 


-y\ 



I y)fj_] rn[\x-y\-\x-z\-\z-y\] 



(4.17) 



p{w) ~ — ==- — — — — exp[m(|x -y\-2^{x- yf/A + w"^)] for m ^ oo. 

y 2t{[[x — yy + w^\-^ 



(4.18) 

Obviously only w of order ^J\x — y\/m are now of relevance in the weight function. 
Therefore one may expand the square roots and gets the simple result 

N 16 , m ,3/2 2 r 2mw^ , . /, 
pH ^(^3^r^ exp[-^-^] form ^00. (4.19) 

This is just the picture which one expects for large mass. Only stream lines near 
the straight line connection essentially contribute, the maximum of p{w) is at Wmax — 



X - y\/{2m). 

If the mass is large enough, such that the variation of the gluon field in transversal 
direction becomes negligible, all line integrals give the same contribution. The weighted 
superposition over the paths can then simply be replaced by the path along the straight 
line connection. This means that one has effectively a three-dimensional (5-function in 
transversal direction. The situation looks now similar to the case of the static propagator 
but with an important difference. While the static propagator (1.8) singles out a special 
reference frame, our propagator is manifestly Lorentz covariant. It is the vector {x — y)ix 
which specifies the direction of propagation. 

This has a simple consequence. In the limit of large mass, the scalar function s in 
(3.9) becomes 



s ~ -\{7;-T^ — r)3/2e-H^-2/IPexp[-i5 r Adz], (4.20) 
2 27r X — v\ Jv 



1 , m 

y\' ~ ~ ' Jy 

with the Euclidean path along the straight line connecting x and y. 

Consider now the term s plus the first two lines of v^^^/j, in (3.11). This sum can be 
written as 



-*E 



E 



—(D^(x)s-sD, {y)h^ 



- ( 



m 



2n\x — y\ 



^3/2; 



1 + 7f (a; - y)f./\x 



y]_^-m\x-y\ 



Pexp[-ic/ r Adz]. (4.21) 

Jy 
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In contrast to the e~™l^''~^"'l of the static approximation which falls off with the euclidean 
time difference, our e""^'^"^' falls off with the euclidean distance. This means that the 
Hamiltonian is correctly given by the full relativistic energy. Furthermore we also get 
the correct projection operator for the 7-matrices. 

We finally discuss the terms a^,t^^ in (3.12), (3.13) which can be treated rather 
simply. We specialize to the case x = y and put X4 — 2/4 = T > 0. In the large mass 
limit we can replace p(w) by 5'^^)(w). The s'-integrations can be written as follows: 

u{z{s\ 0))ds' = u{At = ^dt = - 3 dt. (4.22) 

In the second step we used the 4-component of (3.2), in the third step we introduced the 
asymptotic formulae (4.16), (4.17) for the special case where z lies on the line connecting 
X and y. In we only need to consider the index A = 4 which gives a leading factor 
(— m) from the differentiation of e"'""^. All other contributions are suppressed by higher 
powers of 1/m. In t^^a^i, only spatial indices ^ = = n survive if we concentrate on 
the diagonal part of the Dirac matrices. The transition to 2 x 2-matrices than gives 



in the axial vector (3.12): a^7^7^ =^ emnikFnkl^lm =^ emnkFnkCTm =^ '2BmCrm, 

in the tensor (3.13): t'^'^o-^z. Fmn(^mn emnkFmn<yk 2Bkak. (4.23) 

Both contributions give identical magnetic field insertions which add up. 

The resulting expressions have to be combined with the leading term (4.21) where we 
can drop the projector (1 + 7|')/2. This gives a spin dependent expression of the form 

- (^)3/2e-™^Pexp[-z(7 r Adz](l + ^ / B(t)sdt). (4.24) 
2ttT Jy m J 

If we take the product of the quark- and the antiquark propagator which arises in the 
four-point Green function, focus on the product term of the two magnetic field insertions, 
and extract the Hamiltonion in the usual way we immediately obtain the spin spin 

and the tensor terms with the correct representations for the corresponding potentials 
V4 and V3. In the static approach these terms only arise as higher order corrections. 

Spin orbit terms are momentum dependent and therefore involve moving quarks. 
These are, of course, contained in our formalism, but the derivation is slightly more 
complicated. Spin independent corrections are obtained even harder. Therefore we will 
not discuss those in this first application, but be content with the simple and correct 
derivation of spin spin and tensor terms given above. 



5 Conclusions and outlook 

Let us first compare our representation (2.3), (3.9)-(3.13) with the static approximation 
(1.8). The static approximation is, drastically speaking, completely wrong everywhere. 
It is completely wrong for x 7^ y where it vanishes, but it is also completely wrong for 
X = y because it has a 5-function there which is not present in the exact propagator. 
These features survive if one treats the neglected spatial part i'j"^Dm in the field equation 
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as perturbation. In any finite order of tliis perturbation tlie approximated propagator 
vanislies for x 7^ y, while higlier order derivatives of (5*^^''(x — y) appear. The miracle 
that one can nevertheless obtain useful results from this propagator is due to the fact 
that the perturbation series turns out to become an expansion with respect to 1/m. As 
long as < /m? > is small, the results are reliable. 

The propagator proposed in the present work is manifestly covariant and appears to 
have a reasonable structure everywhere. This advantage is payed by a more complicated 
form, which, however, looks very natural physically. Not only one path ordered integral, 
but a whole set of them contribute. We saw how the paths near the straight line con- 
nection dominate for large mass. We believe that it will also be possible to get useful 
information for finite mass. For an investigation of the quark-antiquark interaction one 
should start, as usual, with the gauge invariant qq four-point Green function and inte- 
grate out the fermion fields. Instead of the familiar Wegner- Wilson loop one will now 
obtain a superposition of loops where the straight paths in time direction are replaced by 
the stream lines making up our propagator. Quite a lot of knowledge how to treat such 
loops has been accumulated by various authors which can be used for this investigation. 

A comparison with the Feynman-Schwinger representation (see e.g. [Q) is also in- 
structive. This representation of the propagator is formally exact and has essentially the 
form of a quantum mechanical Green function. It can therefore be written as a path 
integral. In the literature one also finds approximate path integral representations for 
the propagator and the quark-antiquark kernel valid up to order In both cases one 

has, as usual in this formalism, a sum over all paths, which is, conceptionally as well as 
technically, a rather delicate concept. Even for rather simple situations the path integral 
cannot be evaluated. In our case, on the other hand, we have only line integrals along 
a well defined set of field lines. This is a much simpler situation. Our representation 
stands just between approximations which involve a single path only and those requiring 
a sum over all paths. 

Besides the application of the present propagator for large as well as for finite mass, 
there is another topic which should be worked out. This is the systematic improvement 
of our propagator. We don't have a simple differential equation for it, as it is available in 
the case of the static propagator. Therefore one probably has to improve higher orders 
of perturbation theory directly and transform them into gauge covariant expressions in 
an analogous way as done here for the first order. We emphasize, however, that such 
an improvement does not appear necessary for many purposes. The present form gives 
already the correct relativistic energy of a free particle together with the correct spin 
projectors. Furthermore we have seen that it gives the correct spin-spin and tensor forces 
for heavy quark-antiquark systems. To get these "relativistic corrections" from the static 
propagator one has to make quite complicated manipulations. We expect that all other 
relativistic corrections can also be obtained with some more effort. 

Therefore there are good reasons to believe that the suggested expression for the 
propagator will turn out quite useful already in it's present form. 

Acknowledgement: I thank H. G. Dosch for reading the manuscript. 
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